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q . Abstract 

£>y The role of parallel ion motion for zonal flow generation in ion-temperature- 
^ ! 

, Ph . gradient (ITG) mode turbulence is investigated with focus on the ef- 

, fects of acoustic modes and toroidicity on the zonal flow. One pos- 

£f) ' sible reason for the weak suppression of ITG turbulence by zonal 

^t- : 

■ flows found in experiments in the Columbia Linear Machine (CLM) 

, [Phys. Plasmas 13 055905 (2006)] might be due to the small toroidic- 

o : 

Q\ ' ity (e n = 2L n /R) in the experiment. The zonal flow is often directly 

o : 

dependent on the ITG mode and the coupling of zonal flow to acous- 

•i-H . 

^ ■ tic modes and hence is directly affected by any change of the relevant 

h : 

parameters. The model consists of the continuity, temperature and 
parallel ion momentum equations for the ITG turbulence. The zonal 
flow time evolution is described by a Hasegawa-Mima like equation 
and a fifth order zonal flow dispersion relation is derived. The results 
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are interpreted in terms of quality of zonal flows, i.e., the ratio of 
growth rate and real frequency (Q = Q, IM /Q, RE ). It is found that the 
quality of the zonal flow rapidly decreases with decreasing toroidicity. 
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I Introduction 



In recent experiments in the Columbia Linear Machine (CLM) pQ it has 
been shown that the turbulence suppression of zonal flows is rather weak [2J . 
The suppression of drift wave turbulence by quasi stationary zonal flows 
has been a major topic during the past years in magnetic fusion research. 
The zonal flows are poloidally and toroidally symmetric (kg = 0, k\\ = 0) 
and radially inhomogeneous (k r ^ 0) flow structures in toroidal plasmas [3]- 
[5]. Zonal flows are widely believed to be the main stabilizing mechanism 
for turbulence and one of the main keys to gain access to high confinement 
regimes in present and future magnetic confinement fusion devices such as 
ITER [6] . The generation of zonal flows from drift wave turbulence has been 
extensively studied analytically [3J, [?]- [18] and in computer simulations 
using gyrokinetic [19]- [23] and advanced fluid models [25]- [29] . 

The Ion- Temperature-Gradient (ITG) driven mode is believed to be re- 
sponsible for the anomalously high heat transport in fusion devices and it 
has been shown that it is one of the main ingredients in generation of zonal 
flows [15]- [IB]. The ITG mode has by now been investigated in many dif- 
ferent contexts and situations [30]- [36]. Anisotropic effects in perpendicular 
and parallel ion temperature gradient are out of the scope of the present 
analysis for the zonal flow generation as well as for the linear state [33]. The 
main features of the ITG mode is captured by the more simple model with 
isotropic ion temperature and temperature gradient [37J. Other features of- 
ten found in the CLM experiment, e.g impurities, do not change the ITG 
mode qualitatively, however they may contribute though an impurity driven 
mode [38J that in turn may drive a zonal flow. It has also been show that the 
direct interaction of E x B flows with zonal flows are rather weak [13], [18] . 
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In this particular paper both the slab like and toroidal branch of the ITG 
mode will be considered as the generating mechanism for the zonal flow. 
Here, the ITG mode physics is based on the equations for ion continuity, 
ion temperature and parallel ion motion [21]- [31]. The system is closed by 
considering Boltzmann distributed electrons and using the quasi neutrality 
condition. 

The secondary generated zonal flow time evolution is governed by a 
Hasegawa-Mima like equation. The zonal flow excitation is studied through 
coupling of coherent modes. The drift wave is considered to be a monochro- 
matic pump wave which may couple to the zonal flow through symmetric 
sidebands and this gives a system of four coupled waves [8]- [13] . The corre- 
sponding system may be reduced to a fifth order dispersion relation for the 
zonal flow. It was shown in Ref. [12] that the coherent mode coupling mod- 
els agree well with the corresponding model derived using the wave kinetic 
equation approach. 

To estimate the zonal flow growth rate and real frequency the mode cou- 
pling saturation [39] level is employed. In Ref. [23] simulations of ITG modes 
in CLM were performed and a good agreement between the simulations and 
the experiments were found on the saturation level. This particular finding 
validates the use of the mode coupling saturation level used in this study. 

However, the efforts of comparing analytical estimations with measure- 
ments are rather limited. The effects of coupling between the acoustic modes 
and the zonal flows are mostly unexplored as well. It has been indicated in 
gyrokinetic turbulence simulations that small density gradient scale lengths 
(L n ) increases the transport [22] which encourages the present theoretical 
study. These issues will be discussed in the present study with an emphasis 
on the differences in the generation from the slab like and toroidal ITG mode 
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branch. 

The results are also interpreted in terms of quality of zonal flows i.e. the 
ratio of zonal flow growth rate and real frequency (Q IM /Q RE ) [18] . The 
toroidally driven zonal flows exhibit a much larger ratio of growth rate and 
frequency compared to slab like driven flows (factor of three difference). The 
slab ITG generated zonal flow is only weakly dependent on rji. Moreover, 
the quality of zonal flows is rather insensitive to rji, in contrast to the results 
reported using the Wave Kinetic Equation (WKE) modeling. In many as- 
pects the results found comparing the WKE modeling and the coherent mode 
coupling are qualitatively similar [12]. In the WKE modeling a strongly de- 
creasing zonal flow quality with increasing rji was found [IE], however in 
the present model the decrease of zonal flow quality is visible for larger rji. 
It is also important to note that the difference in quality for slab like and 
toroidally generated zonal flows is approximately a factor of three. It is in- 
dicated that there is a transition from a state with stationary zonal flows to 
a state with oscillatory zonal flows. In all parameter scalings the coherent 
mode coupling model agrees well with a model derived using the WKE model 
for the same background physics. 

The paper is organized as follows. In section II the toroidal and slab ITG 
driven modes are presented and in section III the equations for the zonal flow 
are discussed. Section IV, contains results and discussion and finally there is 
a summary in section V. 
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II Toroidal and slab ion-temperature-gradient 
driven modes 



The model for the ITG driven modes consists of the ion continuity, ion tem- 
perature and the parallel ion momentum equations [21]- [33]. The system 
of equations is closed by using quasi neutrality with Boltzmann distributed 
electrons. The effects of magnetic shear and finite beta are neglected in the 
present study. It has been found earlier that the effect of parallel ion motion 
on the ITG mode growth rate is rather weak [30]. It has recently been found 
that magnetic shear may modify the zonal flow generation significantly for 
the Trapped Electron Modes (TEM), however, the effects of magnetic shear 
on the ITG mode are modest and are also out of the scope of the present pa- 
per [17J. In this section the linear regime of the ITG mode will be discussed 
and the basic equations are given by; 
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Here [A, B] = (dA/dx)(dB/dy) — (dA/dy)(dB/dx) is the Poisson bracket 
with the additional definitions n = Sn/n , 6 = e66/T e , % = 5Ti/T i0 as 
the normalized ion particle density, the electrostatic potential and the ion 
temperature respectively. In the forthcoming equations r = Tj/T e , v+ = 
p s c s y/L n} p s = c s /Vt ci where c s = JTjrrhi, Vt ci = eB/rriiC. We also de- 
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fine Lj = — (dlnf /dr)~ , r\i = L n j Ly., e n = 2L n /R where R is the major 
radius and q.{ = r(l + i]i). The perturbed variables are normalized with 
the additional definitions n = (L n / p s )5n/n , (ft = (L n /p s e)6<p/T e , % = 
[L n j p s )5Ti/T i0 , i>j|| = (L n / p s )5vi\\/c s as the normalized ion particle density, 
the electrostatic potential, the ion temperature and the parallel ion veloc- 
ity respectively. The perpendicular length scale and time are normalized to 
p s and L n /c s , respectively The geometrical quantities are calculated in the 
strong ballooning limit (8 = 0, g± {6 — 0, k) = 1/k where (8) is defined 
by u D (6>) = uj^nQi {8)) [33], [3l]. The effects of safety factor and magnetic 
shear can also be accommodated by the factor gi [15] . The linear dispersion 
relation coming from Eqs 1, 2 and 3 is, 



Here the fields for the linear ITG mode are expanded according to £dw( x , V, t) = 
^ DW e i( ~ k * x + k y v + k w z -^ + c.c. for (n^, 4>,f h Vi\\). The dispersion relation for the 
linear ITG mode including parallel ion motions effects is of 3rd order, however 
if fc| effects are neglected the dispersion relation is transformed into the pre- 
vious result [15J. The effects of parallel ion motion on the real frequency and 
the growth rates are small. The important differences are that the real fre- 
quencies and growth rates are significantly changed for small toroidicity and 
in addition the ^-threshold is decreased. In this particular case it is of inter- 
est to monitor the zonal flow generation from toroidal and slab like ITG mode 
turbulence. In the limit of small e n = 2L n /R and small Finite Larmor Radius 
(FLR) effects an approximate solution may be found, oj 2 slah ~ — fcjjr^j — 2/3), 




(4) 
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in analogy with Ref. [30]. In comparison with the model in Ref. [3D], the 
present model yields a positive threshold value rjuh = 2/3 for the slab like 
ITG mode. In the toroidal ITG mode branch approximate solutions are close 
to those found without parallel ion motion e.g. Eqs 8 - 10 in Ref. [TS] (given 
here for convenience), 



where u = oj r +ij. The origins of the slab and toroidal ITG mode instabilities 
are inherently different: the slab like branch is driven unstable by parallel 
acoustic waves; while the other branch is driven by toroidicity |40j . 

Ill The model for zonal flow generation 

In this section the derivation of the zonal flow dispersion relation will be 
explained briefly. The method used has been covered extensively in previous 
papers [3], [8]- [13J. In this paper the effects on the zonal flow growth rate 
and real frequency of toroidicity and the coupling to the acoustic modes will 
be investigated. There are some limitations to this derivation of the zonal 
flow generation namely the non-linear terms in the equation for tempera- 
ture ([(f), Tj] = 0) and in the parallel ion momentum equation are neglected 
([(f), Vi\\] = 0). The procedure is to start with a monochromatic pump wave 
(defined above) that drives the zonal flow through a modulational coupling 
via the sidebands. The fields for the zonal flow and the sidebands are ex- 
panded as follows, 




CzF(x,y,t) 



(8) 
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i±(x,y,t) = ^ ±e i((l^ x )x±k y y±k llZ -oj ± t) + cc ^ 

The secondary generated zonal flow is defined as £zf (4>,Ti remembering 
that the zonal flow does not have a density component), £± are the possible 
sidebands (n i? 0, Tj). The radial, poloidal and parallel wave numbers are 
(k x ,k y ,k\\), uj is the complex frequency of the drift wave, q x is the radial 
wave number, Q is the complex frequency of the zonal flow. The frequencies 
couple as uj± = Q ± uj. 

The equation for the time evolution for the zonal flow, that is assumed 
to be Hasegawa-Mima type, is expanded using the above fields, 

- ^1<Pzf - nV A x <j>zF = [<P, Vi0] . (10) 

In addition the sidebands are to be determined. The sidebands evolve 
according to Eq. 1 and also identify the coupling between the drift wave 
and the zonal flow as the sideband component. This gives a problem with 
four coupled waves where the zonal flow seed may be unstable and grow 
exponentially. This gives the equation system, 

(itiql- nt)(f>ZF = k y q 2 x [a + <f)* DW (j) + + a_(f) DW (t)-} (11) 
_ k k± uj 4>dw4>zf / 12 s 

0_ = -iPq x k y k ± -u;- — — — — — t^(-13) 

oj_(oj_ + k y - e n gi(3k y ) - kj(3 

In deriving these equations the FLR effects are considered to be small. Here 

a± = 2k x ±q x (14) 
k ± . = kl + k 2 y -q 2 x (15) 
(3 = i + t + t5 (16) 

(17) 



5, 



uj + \rt n gk y 



9 



It should be observed that if k\\ — > 0, Eqs. 11 - 13 are reduced to the 
corresponding equations derived without parallel ion motion [12]. It is also 
important to observe the fact that at q x = the Reynolds stress term vanishes 
and the drive for the secondary generated zonal flow disappears. Using Eqs. 
11 - 13 it is now possible to derive a fifth order dispersion relation for the 
zonal flow component. 

(ft + z/i^)0 + ft_ = [K + {tt + uj)ti^ + K_{tt-uj)ti + ]\<j) DW \ 2 (18) 



= (tt + uj)(Q + a) - kp (19) 

= (ft - - a) - kp (20) 

K+ = k y k ± _/3a + (21) 

K_ = k y k ± _(3a- (22) 

a = uj - k y + e n gik y f3 (23) 



A detailed study of the dispersion relation shows the explicit effect of 
parallel ion motion on the zonal flow. Expanding the dispersion relation 
around the solution ft f° r k\\ =0, where the dispersion relation is reduced 
to third order, with Q = Q + 5Q and solving for 5fl, 

~ W$y {K - " K +) U r\<t>DW?. (25) 

Here it is assumed that the zonal flow growth rate is much larger than the 
real frequency (ftg M >> ^o E )- The important observation from Eqs. 24 
- 25 is the rapid variation with the parallel wave number in the change in 
the zonal flow growth rate and the real frequency in the case of slab like 
ITG mode generated zonal flows. The slab like generated zonal flow growth 
rate and real frequency change in the same manner as the slab ITG mode. 
This also means that the zonal flow is driven by the acoustic waves through 
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the slab like ITG mode instability. The Eqs. 24 - 25 also indicates a weak 
dependency of the parallel wave number on 5fl for the zonal flows generated 
from toroidal ITG mode turbulence. 

The fifth order dispersion relation (Eq. 18) is solved numerically and 
special attention is given to the behavior of zonal flow generation as a function 
of rji, e n and k\\. The results are also contrasted to previous results obtained 
in models excluding parallel ion motion. 



IV Results and discussion 

In the present section the solutions to the zonal flow dispersion relation are 
discussed and compared with previous findings. The zonal flow dispersion 
relation Eq. 18 is solved numerically. In the numerical analysis it is as- 
sumed that the mode coupling saturation level is reached for the drift wave 
turbulence [39] 

7 1 

</W = — j—j— ■ (26) 

Of particular interest in the present work is the effect of toroidicity on the 
zonal flow generation. In this analysis this effect is included in the parameter 
e„ where the slab like ITG mode is found for e n = 0. The results are presented 
both in direct scalings and as the quality of zonal flows. The quality of 
zonal flows is defined as the ratio of the growth rate and the real frequency 

(Q = n IM /n RE ). 

First the effects of e n on zonal flow generation will be discussed. The 
quality of the zonal flow (Figure 1) with the parallel wave number (H) as a 
parameter is investigated. The parameters are r = 1, k x — k y — q x — 0.3, 
H — 0, T)i — 4.0, fcy = 0.1 (asterisk) and fcjj = 0.0 (squares). Shown in Figure 
1 is also the corresponding quality for the ITG mode fcif = 0.1 (plus). The 
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Figure 1: The zonal flow quality (Q IM /Q RE ) as a function of e n with the 
parallel wave number (fcjj) as a parameter is displayed. The parameters are 
r = 1, k x — k y — q x — 0.3, // = 0, rji — 4.0, fc| = 0.1 (asterisk) and 
fcjj = 0.0 (squares). The corresponding quality for the ITG mode is also 
shown k?, = 0.1 (plus). 

parameters used in this paper are chosen as realistic values of a tokamak and 
for easy comparison of the results. The Figure shows a rapidly increasing 
quality of zonal flows with increasing e n = 2L n /R. A large negative value 
of the quality is more beneficial. The qualitative behavior of the model in 
a variation of the parallel wave number is similar except that for small or 
vanishing e n the slab-like mode generated zonal flows are captured. This is 
clearly indicating a stronger suppression effect of toroidally generated zonal 
flows. This result is to be compared with the absolute values of zonal flow 
growth rate and real frequency shown in Figure 2 in Ref . [12J . The quality of 
the ITG mode shows a significantly different qualitative behavior (opposite 
trend) than the corresponding values for the zonal flow. 

Second, the slab and toroidally generated zonal flow qualities as a function 
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of rji are compared with the parallel wave number (fcjj) as a parameter in 
Figure 2. The parameters are r = 1, k x — k y — q x — 0.3, fi = 0; e n = 
1.0 and fcn = 0.0 (asterisk), A;^ = 0.1 (square); e n = 0.0 and fc| = 0.1 
(diamond). The figure shows that the quality of zonal flows is qualitatively 
rather insensitive to rji, in contrast to the results reported in Ref. [IE]. The 
model in Ref. [18] is derived using the Wave Kinetic Equation (WKE) and the 
results are approximately valid in the long wave length limit (k\ « 1 and 
fc| << 1). In many aspects the results found comparing the WKE modeling 
and the coherent mode coupling are qualitatively similar [12], [15]- |18j . In 
the WKE modeling a strongly decreasing zonal flow quality with increasing 
7]i was found, however in the present model the decrease of zonal flow quality 
is visible for larger rji (not shown in the picture). It is also important to 
note that the difference in quality for slab like and toroidally generated zonal 
flows are approximately a factor of three. The results are shown for values 
well above the linear threshold (r] iths i ab = 2/3 and rjuhtor ~ 1-8). 

In Figure 3 the zonal flow growth rate (positive) and real frequency (neg- 
ative) (normalized to the ITG mode growth) as a function of i]i with e n = 0.0 
(diamond) and e n — 1.0 (square) are displayed. The other parameters are 
r = 1, k x = k y = q x = 0.3, /i = and kj = 0.1. The figure shows 
that the zonal flow is exhibiting the same behavior as the linear growth 
rate in scaling with rji. For all values of rji the toroidally generated zonal 
flows have Q IM /jitg > 1 ; whereas the slab like generated zonal flows have 
Q IM /lire ~ 1- The modulus of real frequency (Ifl-^l) of the slab like gen- 
erated zonal flows are, in comparison, much larger than the values from the 
corresponding toroidally generated zonal flow real frequency. The results 
support the findings that turbulence suppression is small in linear devices 
with small curvature (in the present paper represented with small e n ). 
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Figure 2: The zonal flow quality (Vt IM /Q RE ) as a function of rji for slab like 
and toroidally generated zonal flows is shown. The parameters are r = 1, 
k x = k y = q x = 0.3, /i — 0; e n — 1.0 and fcjj = 0.0 (asterisk), fcjj = 0.1 
(square); e„ = 0.0 and fc| = 0.1 (diamond). 

1.5| ' ' ' ' ' ' 

1 0000000000000000000000000000000000 < > 

0.5- 




Figure 3: The zonal flow growth rate (positive) and real frequency (negative) 
(normalized to the ITG mode growth) as a function of rji with e n = 0.0 
(diamond) and e n = 1.0 (square) are displayed. The other parameters are 
r = 1, k x — k y — q x — 0.3, /i — and fc| = 0.1. 
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Figure 4: The zonal flow quality (Q IM /Q RE ) as a function of k\\ with e n = 0.0 
(diamond) and e n = 1.0 (square) is shown. The other parameters are r = 1, 
k x = k y = q x = 0.3, ji = and ^ = 4.0. 

Third, the scaling of the zonal flow quality (Vt IM /Vt RE ) as a function of 
parallel wave number (/c|) is displayed in Figure 4. In the figure it is shown 
that the quality for the slab like generated zonal flows is quite insensitive to 
a change in the parallel wave number. This is due to a similar behavior of 
zonal flow growth rate and frequency with parallel wave number e.g. Figure 5. 
The quality of the toroidally generated zonal flows improves with increasing 
parallel wave number. 

The zonal flow growth rate (positive) and real frequency (negative) (nor- 
malized to c s /L n ) as a function of k\\ with e n = 0.0 (diamond) and e n = 1.0 
(square) are displayed in Figure 5. The other parameters are r = 1, k x = 
k y = q x = 0.3, ji = and r/i = 4.0. A rapidly increasing zonal flow growth 
rate and real frequency is found for small parallel wave numbers and this 
result is corroborated with the approximate Eqs. 24 - 25 of the change in 
zonal flow il. This also shows that the slab like generated zonal flows are 
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Figure 5: The zonal flow growth rate (positive) and real frequency (negative) 
(normalized to c s /L n ) as a function of k\\ with e n = 0.0 (diamond) and 
e n — 1.0 (square) are displayed. The other parameters are r = 1, k x — k y — 
q x = 0.3, fj, = and rji = 4.0. 

driven by the acoustic waves through the ITG mode. This is the reason for 
the comparatively larger real frequency of the slab like driven zonal flow. 
The growth rate and real frequency for the toroidally generated zonal flows 
are only weakly dependent on the parallel wave number. 



V Summary 

In this paper, a fifth order algebraic dispersion relation for zonal flow gen- 
eration is derived including the effects of parallel ion motion. The emphasis 
is on investigating the difference in of zonal flows generated from slab like 
and toroidal Ion- Temperature-Gradient (ITG) mode turbulence. This is of 
special interest since it has been found recently that the suppression of tur- 
bulence by quasi stationary zonal flows are weak in linear experiments while 
it is expected that they significantly contribute to turbulence suppression in 
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toroidal devices. 

The physical model for the ITG mode are the equations for continuity, 
temperature and parallel ion motion. The electrons are considered to be 
Boltzmann distributed and the system of equations is closed by the quasi 
neutrality condition. The secondary generated zonal flow time evolution is 
governed by a Hasegawa-Mima like equation. The zonal flow excitation is 
studied through coupling of coherent modes. The drift wave is considered 
to be a monochromatic pump wave which may couple to the zonal flow 
through symmetric sidebands, giving a system of four coupled waves. The 
corresponding system may be reduced to a fifth order dispersion relation for 
the zonal flow. 

In this study the results are also interpreted in terms of quality of zonal 
flows i.e. the ratio of zonal flow growth rate and real frequency (Q IM /Q RE ). 
The important finding is that the toroidally driven zonal flows exhibit a 
much larger ratio of growth rate and frequency compared to the slab like 
driven (factor of three difference), i.e. they have a higher quality. In most 
cases in scalings with rji the zonal flow growth rate and frequency follow 
the same qualitative behavior as the linear growth rate and real frequency. 
More importantly the quality of zonal flows is therefore rather insensitive 
to a scaling in rji, in contrast to the results reported using the Wave Kinetic 
Equation (WKE) modeling. In many aspects the results found comparing the 
WKE modeling and the coherent mode coupling are qualitatively similar. In 
the WKE modeling a strongly decreasing zonal flow quality with increasing 
rji was found, however in the present model the decrease of zonal flow quality 
is visible for larger rji (not shown explicitly in this paper). In all parameter 
scalings the coherent mode coupling model agrees well with a model derived 
using the WKE model for the same background physics. 
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It is also shown that the slab like generated zonal flows are driven by the 
acoustic waves through the ITG mode. This is the reason for the compara- 
tively larger real frequency of the slab like driven zonal flow. The growth rate 
and real frequency for the toroidally generated zonal flows are only weakly 
dependent on the parallel wave number. 

To this end, in terms of quality, the toroidally driven zonal flows have 
higher values compared to the slab like generated zonal flows such as those 
found in CLM. The slab like generated zonal flows are inherently more os- 
cillating zonal flows and hence the turbulence suppression found is much 
smaller. [2J It should also be noted that the toroidally generated zonal flows 
exhibit large values of the growth rate (fl IM /jitg > 1) compared to the slab 
case. All these results support the findings that the turbulence suppression 
is smaller in linear devices with small curvature (in the present paper repre- 
sented with small e n ). There are still some unknown factors e.g. the effects 
of anisotropic ion temperature and temperature gradient that may shift the 
result. This will be the subject of a future paper. 
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